Hadron production at the LHC: Any indication of new phenomena 
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We confront soft Pomeron and gluon saturation models with the first LHC data on inclusive 
hadron production. We claim that while the first type of models are not able to describe some 
part of the LHC data, the Colour-Glass-Condensate (gluon saturation) approach gives an adequate 
description of the data. Here, we compare our published predictions with the recently available 
7 TeV data. We firmly believe that if further experimental measurements confirm that the gluon 
saturation works, it will be a major discovery. 



I. INTRODUCTION 



As it was expected the first data from the LHC are on soft interaction at high energies such as the total and 
diffractive cross-sections, the hadron inclusive production and so on. Therefore, it is a proper time to review our 
\Q ' understanding of these processes. 

For four decades, the main tool for description of soft interaction has been the approach based on Reggeons and 
Pomerons and their interactions. The first good news is that actually for the first time in these four decades we 
' (— i '[ obtain Pomeron from theory. Today Pomeron is not a plausible assumption as it was in the 70s, it is not the object 
O ■ of successful high energy phenomenology but it comes out naturally from the first theory of strong interaction: N=4 
Super Yang Mills (N=4 SYM) [L|. N=4 SYM together with AdS/CFT correspondence allows us to study the regime 
J^" 1 , of the strong coupling constant Q- For the first time we have a theory which leads to the main ingredients of the 
high-energy phenomenology such as Pomerons and Reggeons, in the limit of strong coupling. On the other hand, 
i—i, N=4 SYM with small coupling leads to normal QCD like physics [U, H| with OPE and linear equations for DIS as well 
as the BFKL equation for the high energy amplitude. 

First, we recall that N=4 SYM has a simple solution for the following set of couplings: 

g s = 4^ = «™ = t4t; r = a ' hxk > 9s « i; but A » 1, (1) 

I/"} . 47T 47rA c 



where R is the radius in AdS§- metric, jV c denotes the number of colours and a' is the slope of the Reggeons 
(a' s» 1 GeV~ 2 ) which is intimately related to the string tension in string theory. The Pomeron which appears in N=4 
SYM [l| at large A, has a intercept and a slope of the trajectory that are equal to 

a F (0) = 2 - JL c4(0) = 0, (2) 

in the limit of < 1. In the next section we will discuss the main property of the Pomeron in N=4 SYM. We will 
introduce available models which incorporate those properties. In the framework of N=4 SYM motivated models for 
soft high-energy interaction, we will then confront these models with the LHC recent data for the inclusive hadron 
production. Our main conclusion is: it is premature to claim that soft model is unable to describe the LHC data, and 
we need to have precise values of the cross-sections of diffraction processes in order to claim so. On the same footing, 
we should improve the Monte Carlo based simulation models in order to include those processes. 

In the last section, we consider the high-density QCD picture which provides an alternative description of soft 
high-energy interaction. We will show that the high-density QCD scenario for inclusive hadron production in proton- 
proton collisions works and gluon saturation reproduces the LHC data for charged-hadron transverse-momentum and 
multiplicity distribution in rapidity and energy [f|. We show that high-density QCD Q predicted 7 TeV data in pp 
collisions [6\. We firmly believe that if further experiments confirm that the gluon saturation works it will be a major 
discovery. There exists some ideas how to simulate the CGC state in nucleus collisions but evidence for the formation 
of the CGC state (gluon saturation) in proton-proton interaction will be a triumph of the high-density QCD. The last 
section is devoted to our predictions for the LHC both for hadron-hadron (pp) and nucleus-nucleus (^4^4) collisions 
[3] • We consider this section as a key part of this manuscript since we believe that comparison with our predictions 
will provide the first hint toward discovery of the new phase of QCD, the CGC at the LHC. 
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II. INCLUSIVE HADRON PRODUCTION IN SOFT POMERON MODELS 



First, let us summarize what we have learned about soft high-energy interaction during the last four decades. 
The phenomenology based on the Pomerons and their interactions has been very successful in description of the 
experimental data. On the other hand all attempts to build a theory of Pomeron and Reggeons and their interactions 
have failed. As it was already mentioned, at present we have two theoretical guides: high energy scattering in N=4 
SYM which gives Pomeron interactions, and matching with perturbative QCD where we can use the BFKL Pomeron 
calculus to obtain the scattering amplitude. Based on these two guides, one can formulate the basic criteria for such 
models: 

(i) Pomeron Ajp ps 0.3 and a' P pa 0. The large intercept of the Pomeron follows from N=4 SYM and its value 
about 0.3 follows from the description of the DIS data |8|. 

(ii) Large Good- Walker component JH since in N=4 SYM the main contribution stems from the processes of elastic 
scattering and diffraction production [if, [l(| ■ 

(hi) Small Pomeron interaction which is of the order of 2/y/X <C 1 in N=4 SYM (Tlj . 
(iv) Only triple Pomeron vertex is essential to p rovide a natural matching with perturbative QCD [iiq . 
At the moment we have only four models [13l - [la | on the market that satisfy these criteria. Unfortunately, the 
inclusive hadron production was predicted in only one model [l3| (see Ref. 
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FIG. 1: Hadron multiplicity from the soft Pomeron model 
divided by <jnd- The experimental data are from Refs. 



at various energy. Dotted curve shows the prediction of Ref. 
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In Fig. [21 we show the typical Mueller diagrams which is used for the calculation of the inclusive hadron production. 
In Pomeron type model, we have two main ingredients: the full Pomeron Green function and &i n defined in Fig. [2] It 
turns out that the main uncertainties stems from the value of <7j n . Various predictions for <7j n can be extracted from 
Table I. 



dN/dy = (1/oin) d a/d y = 1/airi 
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FIG. 2: The Mueller diagram for the inclusive production of hadrons in Pomeron approach. The wave line denotes the full 
Green function of the Pomeron. atot, <J e i, Vsd and add denote the total, elastic, single and double diffraction cross-section, 
respectively. 



The curves in Fig. Q] were calculated using <7j n = <Jnsd- The dotted curve corresponds to Oi n = <Jnd- One can see 
that we need to understand better what kind of inclusive production has been measured. The experimental selection 
of the measured events at the moment apparently depends strongly on the Monte Carlo simulation that has been 



3 





Tevatron (1.8 TeV) 
G(08) G(10) K(07) K(10) 0(C) 


LHC (14 TeV) 
G(08) G(10) K(07) K(10) O(C) 


otat (mb) 


73.29 74.4 74.0 73.9 73.0 


92.1 101 88.0 86.3 114.0 


a el (mb) 


16.3 17.5 16.3 15.1 16.8 


20.9 26.1 20.1 18.1 33.0 


o- sd (mb) 


9.76 8.87 10.9 12.7 9.6 


11.8 10.8 13.3 16.1 11.0 


o-dd(mb) 


5.36 3.53 7.2 13.3 3.93 


6.1 6.5 13.4 12.9 4.83 


ctnsd (mb) 


47.2 46.9 46.8 43.5 47 


61.6 64.1 54.6 51.2 60 


o-jvD(mb) 


41.8 43 41.2 42 


50.1 57.6 41.2 56.2 



TABLE I: Comparison of various soft models: G(08) denotes the GLMM model [l| in which one sums only enhanced diagrams. 
G(10) denotes the model of the Tel Aviv group where a general class of the Pomeron diagrams is summed K(07) and 

K(10) denote two models of the Durham group. The 2007 predictions are from Ref. [3] and the preliminary 2010 predictions 
are taken from the talk of A. Martin at Diffraction 2010 workshop. OS(C) is the model developed in Ref. [la ]. 



v/i TeV 




Pythia6 


Phojet 


GLMM 


0.9 


0-jvd (mb) 


34.4 


40.0 


39.23 


0.9 


cr sd (mb) 


11.7 


10.5 


8.24 


0.9 


cr dd (mb) 


6.4 


3.5 


3.83 


7.0 


ctnd (mb) 


48.5 


61.6 


51.47 


7.0 


o\, d (mb) 


13.7 


10.7 


10.2 


7.0 


(Jdd(mb) 


9.3 


3.9 


6.46 



TABLE II: Comparison of Monte Carlo simulation models with the GLMM model [LI ]. 

used. Table II shows that there is no Monte Carlo code on the market that can consistently describe the diffraction 
production processes. 

Having these two tables in mind, we can conclude that the claim that soft model is not able to describe the LHC 
data is premature and much more work is needed to prove this claim. However, the LHC recent data includes some 
interesting measuments such as the dependence of average transverse momentum of produced hadron on energy and 
multiplicity which cannot be described in framework of Pomeron-type models. 

III. HIGH DENSITY QCD AND HADRON PRODUCTION 

The conclusion from the previous section is that one has to look for a more adequate approach which provides a 
better description of the experimental data at the LHC and will be more closely related to QCD. 

Fortunately we have such an approach on the market: high density QCD [2l[ leads to a completely different picture 
of inclusive hadron production. In this approach, a system of parton (gluons) at high energy forms a new state of 
matter: Colour Glass Condensate (CGC). In the CGC picture, at high energy, the density of partons p p with a 
typical transverse momenta less than Q B reaches a high value, p p oc l/a s 1 (a s is the strong coupling constant). 
The saturation scale Q s is a new momentum scale that increases with energy. At high energies/small Bjorken-x, 
Q s ^> fi where fj, is the scale of soft interaction. Therefore, a s (Q s ) <C 1 and this fact allows us to treat this system 
on solid theoretical basis. On the other hand, even though the strong coupling as becomes small due to the high 
density of partons, saturation effects, the fields interact strongly because of the classical coherence. This leads to 
the a new regime of QCD with non-linear features which cannot be investigated in a more traditional perturbative 
approach. In the framework of the CGC approach the secondary hadrons are originated from the decay of gluon 
mini jets with the transverse momentum equal to the saturation scale Q s (x). The first stage of this process is rather 
under theoretical control and determines the main characteristics of the hadron production, especially as far as energy 
rapidity and transverse momentum dependence are concerned. The jet decay, unfortunately, could be treated mostly 
phenomenologically. 

Actually, such a scenario has passed two critical tests with the experimental data: First, it explains the main 
features of hadron multiplicity in heavy ion- ion collisions at RHIC (KLN papers [.22] ) ; and it gave correct predictions 
for the inclusive hadron production in proton-proton (pp) collisions [H at the LHC at y/s — 7 TeV @. 

The inclusive mini jet cross-section in high-energy pp (or AA) collisions can be calculated within the CGC ap- 
proach via the fc<r factorization [23l - [25j by convolution of two hadrons (or nucleus) unintegrated gluon distributions 
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(f>Q^ A '\xi;kT), depicted in Fig. [3}a, where x\ t 2 = (pt / ' \fs)& ±y Pt and y are the transverse-momentum and rapidity 
of the produced gluon mini jet. The relation between the unintegrated gluon density and the colour dipole-proton (or 
nucleus) forward scattering amplitude N^./^a {xi, rx', b) was obtained in Ref. [24{ which relates the hadron production 
in pp (or AA) collisions to deep inelastic lepton-hadron scattering (DIS) at small Bjorken-rr at HERA. It reads as 
follows 



l) M T ) = M^Wc I ^^^N^ (»,; r T ; b) , (3) 



with notation 



( Xi ; r T - b) = 2N hi{Ai) ( Xl ;r T ; b) ~ (a;*; r T ; b) , (4) 



where rx denotes the dipole transverse size and b is the impact parameter of the scattering. Notice that the relation 
between the unintegrated gluon density and the forward dipole amplitude in the kx factorization is not a simple 
Fourier transformation which is commonly used in literature and also depends on the impact-parameter. The impact- 
parameter dependence in these equations is not trivial and should not be in principle assumed as an over-all factor. 
For the dipole amplitude, we use the b-CGC saturation model [26| which is the generalization of the approach given 
in Refs. [27], [H| and effectively incorporates all known saturation properties [H driven by the small-x non-linear 
evolution equations including the impact-parameter dependence of the dipole amplitude [29]. This model describes 
both the HERA DIS data at small- x [26[, direct-photon production [30j and the inclusive hadron production in pp 
collisions The extension of this model for the case of nuclear target was introduced in Ref. [3] which also give a 
good description of RHIC multiplicity data. 

The main contribution of the kx factorization in the multiplicity comes from px < 2 GeV. For such a kinematic region 
at very low-p-r, we rely on the Local Parton- Hadron Duality principle (3l| . namely we assume that the hadronization 
is a soft process and cannot change the direction of the emitted radiation. This works perfectly in e + e~ annihilation 
into hadrons [3ll l32j | and we believe that this is more preferable than to deal with the fragmentation functions for 
which we have no theoretical justifications at low px- Hence, the form of the rapidity distribution of the mini jet and 
the produced hadron Eq. ([5]) are different only with a numerical factor C, and the transverse momentum of jet and 
the produced hadron are related with a factor (z). In the spirit of the geometrical-scaling property of the scattering 
amplitude, we obtain the charged-particle multiplicity distribution at a fixed centrality but various energies from the 
corresponding mini jet cross-section divided by the average area of interaction a s cx 7r^&| et ^, see Fig. [3]-b. The px 

spectrum of the produced hadron can be then related to the cross-section of the mini jet production in the following 
way: 

^hadron = C " 

Pjet.T 



d2pT - C J *,h[r,]- dj]d2 ^ith PjetjT = PT /(z)\ . (5) 



Notice that kx factorization has infrared divergence. By introducing a new parameter m e f f as mini jet mass which 
mimics the pre-hadronization effect, one can also regularize the cross-section. Therefore, we have only two unknown 
parameters in our model, the overall factor C and the mini jet mass rrijet which are fixed at lower energy. Then our 
results at higher energies and rapidities can be considered as free-parameter predictions. 

Fig. [3] shows our description of the existing experimental data and predictions for higher energies. In this figure, we 
show for the first time the comparison of our prediction |5| with 7 TeV pp data 6] . One can see that the agreement 
is striking. It should be stressed that our approa ch g ives a quite different result for the LHC energies in comparison 
with the Kharzeev-Levin-Nardi (KLN) approach [22| both for pp and AA collisions. The main differences stem from 
the fact that: I) we used a saturation model that describes the HERA data at small-x and has different energy 
dependence and value for the saturation scale. 2) We used a correct relation between the unintegrated gluon-density 
and the forward dipole- nucleon amplitude in the kx factorization, namely Eqs. (|3l4p . 3) In contrast to the KLN, we 
kept explicitly the impact-parameter dependence of the formulation and did not assume that it is trivially factorizable 
as a normalization factor. 4) The relative increase of a s was calculated in our approach while in the KLN approach 
was taken from soft high-energy interactions which is alien to the saturation approach. In both approaches, lower 
energy data was used to fix the overall normalization factor. However, as can be seen from Figs. |3}d, [SJb, we expect 
that the discrepancies between our predictions and the KLN to be more pronounced at higher energies (even more in 
AA collisions). 

In order to obtain the average transverse momentum of charge hadrons, we need also to know the value of the 
average fraction of energy of mini jets carried by the hadrons (z). It is seen from Fig. 0J-a that an average value 
of (z) = 0.48 0.5 is remarkably able to describe the average transverse momentum of charge hadrons in a wide 
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FIG. 3: a) Mini jet production in hadron-hadron collisions in the transverse plane within the kr factorization scheme. The 
impact-parameter between two hadrons is B. b) shows the average impact parameter of the produced mini jet (fej ei ) as a 
function of energy within two rapidity bins, c) The comparison with the experimental data and prediction for dN c h/dy. The 
curves are normalized by data at \fs = 546 GeV p\ . d) Energy dependence of the charged hadrons multiplicity in the central 
region of rapidity r\ = in pp collisions. The theoretical curve (Saturation model LR) is our prediction coming from the 
saturation model [5j. The band indicates about 2% theoretical error. The total theoretical uncertainties is less 6% at high 
energies. We also show the KLN prediction [22l ] with the same error band as ours. Notice that in c panel we have taken a fixed 
mini jet mass m,j e t = 0.4 GeV for all energies while in d panel uncertainties due to the assumption of a fixed energy-independent 
mini jet mas was included in the band. The experimental data are from Refs. [tj [18I - I20I , |33|] . The experimental error bars 
indicate systematic uncertainties. 



range of energies. In order to further test the validity of the value (z) for the mini jets, we show in Figs. @}b, 13} c our 
predictions for the differential yield of charged hadrons in the range \r)\ < 2.4 and at various J 77 1 bins for t/s — 2.36, 7 
and 14 TeV. The experimental data are recently reported from the CMS collaboration 0, It is seen that our 
predictions is in quite good agreement with experimental data at 7 TeV. We recall again that the pre-factor in Eq. §5§ 
is the same as what we already fixed with experimental multiplicity data at low-energy y/s — 546 GeV at 77 = 
in Fig. [SJ-c. Therefore, we have no free parameters in obtaining the theoretical curves in Fig. 2) The fact that our 
model reasonably works at low pr is due to the fact that the saturation scale is rather large at low pr- In our 
formulation, we predicted that the differential yield of charged hadrons has a peak at low pr- The position of the 
peak is approximately at pt — m,j et (z) The experimental data at 7 TeV shown in Figs.HJd indeed confirmed this 
prediction. In the CGC scenario, the gluon saturation scale is proportional to the density of partons. The parton 
density is proportional to the multiplicity and, therefore, one can relate the saturation momentum in the event with 
the multiplicity of the hadrons n. In Fig. [4jd, we show the average transverse momentum of charged hadrons as a 
function of the number of charged particles for events within the kinematic range pr > 500 MeV at various energies. 
The experimental data are from ATLAS for yfs = 0.9 TeV [34|. Our prediction seems also to be in a good agreement 
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FIG. 4: a) The energy dependence of the average transverse momentum of charged hadrons. b,c)The differential yield of charged 
hadrons. The LHC experimental data are from the CMS collaboration H@. d) The average transverse momentum of charged 
hadrons as a function of the number of charged particles for events with n c h > 1 within the kinematic range pr > 500 MeV. 
The experimental data are from ATLAS for -^/s = 0.9 TeV and \rj\ < 2.5 34]. The theoretical curves was obtained for = 
and with the same kinematic constraint pr > 500 MeV at various energies for two value of (z) = 0.48, 0.5 corresponding to the 
dashed and the solid lines, respectively. The mini jet mass is taken m je t = 0.4 GeV in all plots. The normalization is the same 
as in Fig. [3] 



with preliminary 7 TeV data from the ATLAS collaboration (not shown in the figure). 

Finally, in Fig. [3] we show our prediction for Pb+Pb collisions at the LHC [7| . Notice that again similar to the 
case of pp collisions, we have here only two free parameters, normalization factor C and the mini jet mass which are 
fixed at RHIC energy y/s = 200 GeV for — 6% centrality. Therefore, at lower/higher energies than y/s = 200 GeV 
(for various centrality /rapidities) we do not have any free parameter. In Fig. [5}a, we show our predictions at lower 
RHIC energies -y/s = 19.6 and 130 GeV in Au-Au collisions, and also for the LHC energies yfs — 2.75 and 5.5 TeV 
in Pb-Pb collisions for — 6% centrality bin. In Fig. [5}b, we show the energy dependence of dN pp /dr), dN^A/dr] and 
(2/N par )dNAA/dr] at midrapidity rj = for central collisions (where N par denotes the number of participant for a 
given centrality). 

To conclude, we showed that the CGC gives very good descriptions of the first data from the LHC for the inclusive 
charged-hadron production in proton-proton collisions, the deep inelastic scattering at HERA at small Bjorken-a:, and 
the hadron multiplicities in nucleus- nucleus collisions at RHIC. We believe that our predictions for nucleus-nucleus 
collisions at the LHC will be a crucial test of the CGC approach. 
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FIG. 5: a) Pseudo-rapidity distribution of charged particles produced in Au-Au and Pb-Pb central — 6% collisions at RHIC 
and the LHC energies, b) Energy dependence of the charged hadrons multiplicity at midrapidity r\ — in central collisions 
in pp and AA collisions. The theoretical curve Saturation model (LR) is our prediction. The band indicates less than 3% 
theoretical errors. The total theoretical uncertainties is less than 7%. The experimental data are from [a. I18H20, 33, 35]. The 
plots are taken from Ref. 0|. 
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